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Abstract 

We study the phase diagram of scalar field theory on a three dimensional Eu- 
clidean spacetime whose spatial component is a fuzzy sphere. The corresponding 
model is an ordinary one-dimensional matrix model deformed by terms involving 
fixed external matrices. These terms can be approximated by multitrace expres- 
sions using a group theoretical method developed recently. The resulting matrix 
model is accessible to the standard techniques of matrix quantum mechanics. 



1. Introduction and results 

Fuzzy spaces, such as the fuzzy sphere Sp [T], provide an interesting way of regularizing 
quantum field theories |21 E] : These spaces come with an algebra of functions which is finite 
dimensional. Correspondingly, the functional integral appearing in the partition function of a 
quantum field theory on a fuzzy space reduces to a finite dimensional integral. An advantage 
of this approach over a lattice regularization is that it preserves symmetries: the isometrics 
of the classical manifold still have a well-defined action on the corresponding fuzzy space. 

Interestingly, (scalar) quantum field theories on a fuzzy spher^ correspond to hermitian 
matrix models with additional couplings to fixed external matrices [H O [6] . These matrices 
originate from the kinetic term in the action, and they yield an obstruction to applying 
the usual matrix model technology for reducing the partition function to an integral over 
eigenvalues. Nevertheless, one can use group theoretical methods to perform this reduction 
after Taylor-expanding the exponential of the kinetic term in the partition function [5l [6] . 
The rewritten partition function can then be used to derive analytically the phase structure 
of scalar field theory on the fuzzy sphere and to compare the result to the numerical studies 
of [3, El El Uni mi [I2] : Qualitatively, the phase structures match and the position of a triple 
point appearing in the phase diagram roughly agreed in both the numerical and analytical 
results. 

In this paper, we continue the work of O [6] by studying scalar field theories on the 
three-dimensional spacetime IR, x Sp. The resulting model is a special variant of matrix 
quantum mechanics with additional couplings to fixed external matrices. We can use again 
the group theoretical techniques developed in [5]. That is, we Taylor-expand the exponential 
of the spatial part of the kinetic term in the partition function and rewrite the resulting 
series order by order in terms of multitrace expressions. Re-exponentiating these yields the 
partition function of a matrix quantum mechanical model involving the multitrace terms. In 
fact, we can trivially translate the results obtained in [6j to the matrix quantum mechanics 
model. 

Similar to a pure hermitian matrix model, matrix quantum mechanics can be solved in 
the large limit. One approach is to map the system to a non- interacting Fermi gas [13], an 
alternative way is the collective field theory formalism developed in |14j. These techniques 
can still be applied if the potential contains multitrace terms of the matrix field, as done e.g. 
in |15|, [161 [T7] in the 2d gravity context. The model that we obtain will be more general in 
the deformations than the models described in these papers, but it will still be accessible to 
these techniques. 

The purpose of this paper is to complement the numerical results on the phase diagram of 
scalar (/)^-theory on IR, x S"!. found in [HI [l9l [20] by an analytical study. Using the techniques 
mentioned above, we are able to derive many properties of this phase diagram. In particular, 
we confirm the existence of a third phase as compared to pure matrix quantum mechanics. 
Furthermore, we find expressions for the various phase boundaries, and in an important 
region of the parameter space, we can write down explicit formulas for the free energy. In 
the remaining parameter space, we can calculate the free energy for each point using two 
simple numerical operations. 

'^or, more generally, fuzzy projective algebraic varieties 
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This paper is structured as follows. In section 2, we briefly review scalar quantum field 
theories on fuzzy complex projective spaces and show how to turn these into matrix quantum 
mechanics. An expression for the free energy of such a theory on a fuzzy sphere is derived in 
section 3. In section 4, we discuss general features of the phase diagram of this theory and 
quantitative results are presented in section 5. 



2. Fuzzy scalar field theory as multitrace matrix quantum mechanics 

2. 1 . Fuzzy complex projective spaces 

We will work with Berezin-quantized complex projective space as described e.g. in [6] and 
we briefly summarize the associated notions in the following. Consider the space C"'"''^ with 
complex coordinates 0, ...,n. Using the fibration U(l) S^^+i ^ CP", we conclude 

that the set of functions on C""*"^ of the form 

f{z) = ^ jai-afffi-ft ^ yai...Q(,/3i.../3(, ^ q-, (^2,\) 

descends to a subset of the smooth functions on CP" C 5*^"^^ C C""*"^. The Hilbert space 
for Berezin-quantized CP" is the ^-particle Hilbert space in the Fock space of n + 1 
harmonic oscillators. For functions in of the form ()2.ip . the quantization prescription is 
explicitly given by 

/(^) ^ / = E r"-°^''"-^^^at,,-«L|0)(0|a/3,...a;3, • (2.2) 

Note that real functions / are mapped to hermitian operators / G End (=^) and the constant 
function f{z) = c, c G R, is mapped to c • 1 G End (J^). Furthermore, N^^i, := dim(^) = 
'^'^^p ' and a real function / G can be interpreted as an Nn/ x A'^^^-dimensional, hermitian 
matrix. 

The Laplace operator A on S can be lifted to End {M'l) and the lift is given by the 
quadratic Casimir C2 of SU(n + 1) in the representation formed by End (J^). One can show 
that [m 

A/ = C2/ . (2.3) 
Furthermore, we can define an integral operation on End (=^) by taking the trace: 

..^/.Z^IIP (2.4) 

where d^ps is the Liouville measure induced by the Fubini-Study metric on CP". For a 
more detailed exposition of these relations, see e.g. |21j . 



2.2. Partition function of the model 

The action of scalar field theory with quartic potential on R x CPJ- is given by 

= /3 Jdttr (i$(t) {C2 - df) $(t) + r ^\t) + g^\t)) . (2.5) 
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The scalar fields are represented by the time-dependent, A^n/ x ^-dimensional hermitian 
matrices <I>(t) E C°°(1R, End (J^i)) and, as mentioned above, the spatial part of the Laplace 
operator A is given by the quadratic Casimir C2 of SU(n -|- 1). Note that we are working 
with Euclidean time, which implies a different sign in front of the potential compared to 
Minkowski signature. 

Explicitly, we have C2^{t) := [L'\ [L*, <!>]], where the U form generators of SU(n -|- 1) 
acting on .y^. We follow the conventions of [6] and normalize the generators according to 

[Li,Lj]=:i'^fijkLk, tr(Li) = 0, ^L^ = cl1 and tr {LiLj) = ^^2'^, 1 ■ 
k i {n+l) i 

This yields positive eigenvalues for C2, and thus the action (|2.5|) is bounded from below for 
5>0. 

The model (|2.5p defines a one-dimensional field theory of matrix- valued scalar fields. 
Such theories are known as matrix quantum mechanics in the literature. The corresponding 
partition function reads 

Z = j exp(-/35[$]) , (2.6) 

where &^{t) is the measure over functions on IR, taking values in the Lie algebra of hermitian 
matrices of size Nn/ x Nn/. For each value of t, the measure corresponds to the Dyson 
measur^ on the space of hermitian matrices. As in the case of the pure matrix model, we 
can therefore split the partition function into an eigenvalue part ^A(t) and an angular part 

Z = J &A{t)^n{t) exp{-(3S[<^>]) . (2.7) 

One now usually exploits the fact that expressions in the action which consist exclusively of 
traces of polynomials in the $ are independent of the angular part: Diagonalizing <I> according 
to A := ft^ ^Q, where A is the diagonal matrix of eigenvalues of ^ and ft is some unitary 
matrix, yields the simplification 

tr ($") = tr {{nm^)'') = tr (A") . (2.8) 

The fixed external matrices Li originating from the quadratic Casimir, however, present an 
obstacle to rewriting the action in terms of A. In [5] the same problem was analyzed for scalar 
field theory on fuzzy spaces, that is, the dimensional reduction of the model (12. 5p . There it 
was suggested to perform a Taylor expansion of the exponential of the kinetic term in the 
partition function. Order by order, the terms in this expansion can be evaluated analytically 
using group theoretical methods and, when recombined, yield the partition function of a 
matrix model with an action containing multitrace terms. We can apply the same method 
here and directly take over the results obtained in [6j: The Lagrangian of our model is given 

by 

£ = tr ( - ^d^^ + r«>2 + 5^.4) + £casimir , (2.9) 



^i.e. the standard translation-invariant measure induced by the bi-invariant Haar measure on U(A'^„,, 
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where >Ccasimir = tr {^{t)C2^{t)). For the fuzzy sphere GPp, the latter term can be approx- 
imated in the Hmit of large matrix sizes by 

^Casimir = ("^ tr - ^ tr {^f^^ tr ($) + (^1 - A tr ($2) + ^ tr tr ($2) , 

(2.10) 

where we abbreviated N = Ni^i. We thus arrive at a standard one-dimensional matrix 
model, deformed by the multitrace terms contained in £casimir- The corresponding actions 
for CP^ with n > 1 are easily deduced from the results in |Bj, too: They merely correspond 
to replacing the coefficients in (j2.10p by more complicated expressions in N. Since our goal 
is primarily to reproduce the numerical results of [18l [19] , we will focus our attention on the 
case of the fuzzy sphere in the following. 

3. Evaluation of the free energy 

In this section, we use collective field theory [H] to compute the free energy of our one- 
dimensional matrix model. To leading order in the matrix size, this method is equivalent to 
the semi-classical approximation of the model as presented in [13], see also [22] for a nice 
review. 

3.1. The free energy of multitrace matrix quantum mechanics 

We start from the general one-dimensional matrix model with Lagrangian 

/: = i tr ($2) + U2 tr ($2) + U4 tr ($^) 

+ vi2 tr + VM tr ($)^ + V22 tr {^^ f + wu tr {<^ f tr ($2) (3.1) 
= itr($2) + F . 

For simplicity, we decompose the matrix $ into = ^""Ta, where Ta are hermitian 
generators of u(A^), normalized according to tr (raT^) = 5ab- The Hamiltonian corresponding 
to (|3.ip can be written as 

^ = -i tr ($2) + tr («>4) 

(3.2) 

+V12 tr {<^>f + vu tr ($)^ + 1^22 tr {<^^f + wu tr {<^f tr {<^^) . 
We now switch to the collective fields (j){X) G Cc(lR) via the Fourier transform 

m = I ^e^'Hrie--^^'^), (3.3) 

where we inserted a factor of compared to [14j to facilitate the large N limit. Here, the 
collective field </>(A) will turn out to play a similar role as the eigenvalue density in the case 
of the hermitian matrix model. We correspondingly introduce the various moments of 0(A): 

Ck := J dA(/.(A)A'= . (3.4) 
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Note that (|3.3p implies that the collective field (/'(A) satisfies the normalization condition 

co = J dA,^(A) = ltr(l) = l . (3.5) 

The potential V, i.e. the non-derivative terms in the Lagrangian (j3.1|) . is easily seen to 
transform into 

V = ^J dA</.(A)((t;i2Ci + f 14Ci)A + (U2 + V22C2 + tL'l2C?)A^ + U^X'^^ . (3.6) 

The transformation of the kinetic term is technically more involved, see [131. Here, let us 
just note that one eventually arrives at the transformed Hamiltonian Ti = T + V where the 
effective potential V reads as 

V = I dXm{^G\X,cP))+V (3.7) 

with the resolvent G{X, (p) = J d( To determine the collective field (/'o(A) for the ground 
state of the system, we have to minimize the functional 

E{HF, VF, KF,<p) =V + flF (^^ - j dA 0(A)^ 

+ UF (^ci - y dA (A(A)A^ + KF (c2 - y dA </>(A)A 
The variations with respect to (p, fiF, t^F, K'F, ci and C2 yield the system of equations 

= Hf + i^pX + kfX^ - {(yi2Ci + vuc\)X + {u2 + U22C2 + wi2c\)X^ + u^X'^ 
1 = / dA0(A) , ci = / dA</.(A)A , C2= I dX(p{X)X^ , 



(3.9) 



(3.8) 



UF = -V12C1 - Svucf - 2wi2CiC2 , KF = -V22C2 , 

which is solved by standard methods, cf. [HI [23]. We arrive at the collective field 

0o(A) = ( lVnf^F-a,X-a2X^-a,X^) for A G X , ^^^^^^^ 
[ else , 

where 

ai = 2f 12C1 + 4f 14C1 + 2wi2CiC2 , 02 = M2 + 2f22C2 + wucf , 04 = M4 , (3.10b) 

and X C R can be any union of closed intervals such that (j)o{X) = at the boundaries. 
Equations ()3.9p imply that 

J dXM>^)G\X,M = l, j dX4>oW{2{iiF-aiX-a2X^ -aiX'')) = ^ J dX^ , (3.11) 
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and the potential V reads as 




/ 



dA0o(A) 




(/)o(A) + fJ,F - C1W12A - C2W22A^ - 3clvuX + 2C1C2W12X 



) 



Together with (j3.10|) , we arrive at the fohowing expression for the energy (|3.8|) in the ground 
state: 



While all the integrals appearing above can be expressed in terms of elliptic functions as done 
in appendix A, it is not possible to use them to find analytic expressions for /ij? in general. 

3.2. Comments on the semiclassical approximation 

The semi-classical approximation employed in [13] leads to the same result as the collective 
field theory method. Here, one considers the quantum mechanical problem of finding the 
eigenvalues of the A'^-particle Hamiltonian corresponding to the action (j3.1l) . where $ = 
$^ G Mat(c(A''). As in the case of the zero-dimensional matrix model, one can switch to an 
eigenvalue formulation of the model. The Vandermonde determinant arising from this can 
be absorbed in the quantum mechanical wavefunction, rendering it totally antisymmetric. 
The ground state energy of this system is then found by using the standard description in 
terms of free fermions, where the Lagrange multiplier of the collective field theory method 
becomes the Fermi energy. 

The only new ingredient here is the treatment of the multitrace terms: By linearizing e.g. 
tr (<I>^) around the vacuum expectation value C2 = (tr (<I>^)), we find the relation 



cf. [16]. This induces a constant shift of the free energy proportional to (tr ($^))^ as well 
as a doubling of the naive contribution of the term corresponding to tr ($^)^. In |16j . this 
was justified from physics principles. In the previous section, we saw that precisely the same 
modifications arise in our model from integrating out the Lagrange multiplier Kp. 

3.3. The large N limit 

As usual for matrix model^, the large N limit is not unique, but requires fixing of the 
scaling behavior of all coupling constants and the eigenvalues. While mathematically all 
scalings which leave the potential positive definite are equally valid, the resulting models are 
physically very different. One physical constraint one usually imposes is that in the limit 
— )• 00, the coupling constants approach critical values. 

Our model is supposed to provide an approximation for scalar field theory on x S'^, 
and we expect all terms of the potential as well as the kinetic term to survive in the large N 
limit without becoming dominant. This implies that the leading order of all the contributions 
to the action should be of order 0{N^), subleading contributions of order 0{N^^) provide 

similar problem is encountered in quantum field theory on the lattice. 





(3.12) 




(3.13) 
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corrections to the large limit. An overall scaling of the action is irrelevant, as we can 
define the free energy with a corresponding power of N. 

To obtain a homogeneous scaling between the approximations due to the kinetic term 
and the potential, we can use the scaling found in [6j: 

(3-^N'^(3, X-^N'^X, r-^N^r, g-^N^g, (3.14) 



where r and g are the actual parameters of our model (|2.9p . In addition, we have to scale 
t — )■ jj^t to include the right scaling of the temporal part of the kinetic term. Altogether, our 
model corresponds to the general model (|3.ip with the following parameters: 

U2 = l + r, U4 = g, vi2 = -l, ^1"^ ^ ~f ' ^22 = -^, wi2 = ^/3, (3.15) 
and for X £ I, the corresponding collective field reads as 

M^) = (^fiF - (-2ci - pel + ^/3c^C2^ A - (^1 + r - ^/3c2 + '^/^cj^ X^ - gX^ 
We furthermore have the relation 

ai = 2ci(r — 02) . (3.16) 



4. The phase diagram 
4.I. The three phases 

We will exclusively deal with the situation in which tq is negative enough for the potential to 
have a local maximum. Note that in matrix form, our potential is symmetric under 
In the collective field reformulation, we have a symmetry of the potential V under A — )• —A. 
Naively, one therefore expects that symmetric phases are the dominant ones. However, our 
experience with the pure matrix model [6j suggests that we should also allow for a third 
phase. This is further motivated by the numerical findings of \18 \ 119 1 [20]. which also observe 
(at least) three different phases. Depending on fip, we can distinguish three phases: 

I. The disordered phase or single-cut case. Here, fip > and the filling of the eigenvalues 
covers the local maximum completely. There is a single interval X = [— Ai,Ai] over 
which one has to integrate A, and it is symmetric around the origin. This implies that 
the first moment of the collective field vanishes: ci = 0. 

II. The non-uniform ordered phase or symmetric double-cut case. In this case, fip < 
and the Fermi sea of eigenvalues splits up into two symmetric, disjoint pieces. That is, 
there are two intervals I = [— A2,— Ai] U [Ai,A2] with < Ai < A2 as support for the 
integral, which are again symmetric around the origin. We again have ci = 0. 

HI. The uniform (ordered) phase or asymmetric double-cut case. As in phase II, fip < 
and the Fermi sea is split, but this time into two asymmetric pieces. The length of the 
two intervals Xi UX2 = X is different. Note that true uniform ordering is only achieved 
in the totally asymmetric case, in which one of the intervals has shrunk to zero size. 
Here, ci / 0. 
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The actual phase is determined from existence conditions and the fact that the physical 
system adopts the ground state with the lowest free energy. 




Figure 1: The three phases. The solid line describes the potential ^(A) felt by the eigenvalues. 
The dashed lines mark the Fermi energy fip. In phase III, we chose a totally asymmetric 
filling with Xi = 0. The asymmetry of the potential in this phase is due to ci 7^ 0. 

Note that our model can always be treated as ordinary matrix quantum mechanics, where 
the coefficients of the potential are dependent on the moments ci and C2- That is, to study 
our model, we can solve matrix quantum mechanics for a general potential and then derive 
self-consistency conditions on the moments. For this reason, it is possible to evaluate the 
exact location of the phase transition between phases I and II as well as the existence domain 
of phase III in our model analytically. 



4-2. The phase transition I to II 

The phase transition between the single-cut phase and the double-cut phase obviously occurs 
at ^i? = 0. In this case, the elliptic integrals can be performed explicitly. The interval I on 
which <^o(-^) is supported is given by 



-U2-2V22C2 / —U2-2V22C2 
U4 ' Y U4 



and the normalization condition yields 

\/8{-U2 - 2V22C2) 



37ru4 

The second moment C2 is determined by the condition 



V32 {-U2^2v22C2)^ _ (Svr) 



'4 5m| 



(4.1) 



1 . (4.2) 



C2 = 77—2 = ' (4-3) 



where we made use of (j4.2p . Plugging this into (j4.2p yields a phase transition at 

U2 = ,(3^)^(5^4 + 4.22) _ (4 ,) 

10(U4)3 

If U2 is larger than the right-hand side, we are in the single-cut phase, while for U2 smaller 
than the right-hand side, the system is in the double-cut phase. 
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4-3. Existence condition for phase II 

Similarly to the pure matrix model discussed in [6], there is a region of the parameter space 
which is not covered by phase I and where the self-consistency relation for C2 cannot hold true 
in phase II. To see this, we compute C2 for the double cut solution of pure matrix quantum 
mechanics with vi2 = vm = V22 = W12 = for arbitrary parameters 02 < and 04 > 0. We 
then translate to our model by identifying 

2 

02 = 1 + r - -/3c2(a2, a4,//_p(a2, 04)) and 04 = ^4 = 5, (4.5) 

where r and g are the actual parameters of our model. The region of the parameter 
space in which we expect the self-consistency relation to be problematic corresponds to 
< a4/(— 02) <C 1. In this region, the wells in the potential are sufficiently deep to be 
approximated by a parabola. That is, we approximate e.g. the right well as follows 

2{^p-a2X^-a4\')^2fiF + ^ + 4a2(^X-^l^^ . (4.6) 

Both the integrals yielding the normalization condition for (po and the second moment C2 can 
be computed and read as 



/ 



dX (/>o(A)A 



2 _ (ai + 4g^^)(17ai+4g/i^) 
128g2(-a2)3/2 

We can use these results together with (|4.5p to determine the function r(a2): 

r.-l^^-laJ^-,) . (48) 



12y/-a2 3 \g 

Note that the larger \a2\, the better this approximation becomes. In particular, we see that 
for g < ^(3, r grows as 02 becomes more negative. This implies that phase II exists for 
arbitrarily large values of |r| only if 5 > 



4-4- Existence condition for a totally asymmetric phase III 

To simplify our analysis, we will identify the third phase with a totally asymmetric filling. 
That is, I is again a single interval, filling only one of the two wells in the potential. It is 
obvious that the transition from phase I to the totally asymmetric filling has to go smoothly 
through all possible asymmetric fillings, as close to the boundary between phases I and II, 
no totally asymmetric solution will exist. 

The existence boundary for the totally asymmetric solution is another line in our phase 
diagram which can be determined analytically if we neglect the contribution of the odd 
moment ci. As we will see later in our numerical studies, the odd moments decrease the 
depth of the filled well. The bound obtained by this approximation for the existence of a 
totally asymmetric phase III is therefore an upper bound. 

In our approximation, this line is determined by the fact that one of the two wells of the 
potential is filled up to the local maximum and that a further increase in Fermi energy iip 
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would lead to a spilling of eigenvalues into the other well. Consequently, we put = and 
consider the interval 

-U2—2V22C2 



o.A/ '"^;r" ) • (4-9) 



The computation of the existence boundary then proceeds exactly as the computation in the 
previous section. We obtain 



V2{-U2 - 2V22C2)^ , 1 \/8 (-n2 - 2U22C2)2 (Gtt) 

= 1 and C2 = 9 = 1 

3vru4 157ru| 

and the existence domain of the totally asymmetric phase is given by 



(4.10) 



U2 < _(3^)M5n4 + 4.22) _ (4^^^) 
5(2n4)3 

4-5. Pure matrix quantum mechanics 

Before discussing the phase diagram of our model, let us briefly discuss the slightly simpler 
case of pure matrix quantum mechanics, for which vu = f 14 = V22 = = 0. Here, the 
phase transition between I and II occurs at 

= (4.12) 

and the existence boundary for the totally asymmetric cut is 




U2 -5 -4 -3 -2 -1 



Figure 2: The phase diagram for pure matrix quantum mechanics. The solid line describes 
the phase transition and the dashed line the boundary of the existence domain of the totally 
asymmetric cut. The roman numerals describe the actual phases of the system. 
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One can now show by a simple physical argument that the lowest energy configuration is 
always the symmetric double-cut solution, and thus phase III does not exist in pure matrix 
quantum mechanics: In the time- independent ground state and for large N, one would 
expect that moving one eigenvalue from one well to another does not cost or yield energy. 
This automatically implies that the Fermi energy in both wells is the same. To be rigorous, 
one can introduce filling fractions pi and p2 for the left and right wells of the potential 
with pi + p2 = 1 and Jj-_ 4>oW = Pi- One can then determine the true minimum of the 
free energy numerically, which confirms the physical argument. The phase diagram of pure 
matrix quantum mechanics is depicted in figure [2l 

5. Results 

In the following, we put /3 = ^ and perform an analysis of the phase diagram relying on solving 
the system of equations (|3.1U|) and determining the free energy. As we already computed 
the location of the phase transition between phases I and II, it remains to compare the free 
energy on the overlap of the existence domains of phases II and III. 

5.1. Phases I and II 

Recall that the filling of the wells with eigenvalues is symmetric and therefore the odd mo- 
ments of the collective field, and in particular ci, vanish. Equations (|3.10p cannot be solved 
analytically, as they involve elliptic integrals. We therefore apply the following algorithm: 

1. pick a value for 04 = g and 02 

2. evaluate fip for this pair by numerically minimizing |1 — Jj d\ (j)o{\)\ 

3. numerically evaluate the second moment C2 

4. deduce the value of r from 02 and C2 

5. evaluate the free energy at the point {r,g) in the parameter space 

By applying this algorithm to a range of values for 02 and 04, we can perform a sweep of 
the r-5-parameter space. The only restriction here is the existence boundary for phase II for 



The region of the r-(7-parameter space, in which we expect an overlap between phases 
II and III is characterized by large |r| and small g. The potential in this region has broad 
and deep wells, and correspondingly pp is small compared to the depth of the wells. This 
implies that an approximation of the potential by a parabola, like the one used in section 
14. 3. t should work reasonably well. 

Explicitly, we approximate the collective field by 



The zeros of this collective field (/>o(A) defining its support I = [—\r, —Xl] U [Xl, Xr] are 





else . 



for A G X , 



(5.1) 




and 




(5.2) 
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Together with the integrals given in appendix B, we can now compute 



1 /3 2 

po = , C2 = + A^i„ , a2 = 1 + r - - -/3A^i„ . (5.3) 
4-y/a ova 3 



Prom these relations, we derive 



4 3 ' ^ 8A^i„ 

4 48a 8 6^ 3 



(5.4) 



Rewriting the free energy NE in terms of r and g yields a lengthy but analytic expression, 
which is plotted in figure [3l Note that one can easily check the accuracy of the approximation 
by verifying the exact normalization condition as well as the self-consistency condition for C2 
using the approximate values oi jip and C2- This plot is confirmed by the numerical results 
of the first algorithm. 




Figure 3: The free energy in phase II as obtained from the approximation of the potential 
wells by parabolas. 



5.2. Phase III 

Phase III is more intricate to deal with due to the additional appearance of the terms involving 
the odd moment ci. Prom the expression of the free energy ()3.12p . it is clear that the free 
energy in phase III is much larger than the free energy in phase II if C2 and the value of 
the integral are assumed to be roughly the same. We therefore do not expect phase III to 
compete with phase II, but to exist only where phase II does not. This is precisely the case 
in the region of the r- (/-parameter space for which the approximation of the wells of the 
potential by parabolas, and thus (|5.ip . works well. 

For a totally asymmetric filling, i.e. a collective field (/'o(A) with support X = [Al, A/j] we 
obtain 

Pq = , ci = Amin , C2 = —= + A^jj, , a2 = l + r - —= . (5.5) 

z-v/ct 3vo 
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This leads to the fohowing expressions: 



1 „ 2/3 \ a -2 /3 

— 6-3a-^ , g = + 



- 12"'^' + 2/3 + 6aA^i, - Sa^A^j, + 30V^/3A^i, + UafiXj,^ 

24a 

A plot of the free energy in phase III according to this approximation is given in figure SI 




Figure 4: The free energy in phase III as obtained from the approximation of the potential 
wells by parabolas. 



5.3. Discussion 

In the region of the parameter space in which phases II and III coexist, the free energy is 
negative in phase II, while it is positive in phase III. Note that this observation is due to 
our restriction of phase III to a totally antisymmetric eigenvalue filling. Actually, one would 
expect the system to adopt intermediate phases so that the free energy changes continuously 
between the extrema of a symmetric and a totally asymmetric filling. 

Prom the different signs of the free energy in phases II and III it follows that phase 
III is only adopted by the system if phase II is not available. For large |r|, we therefore 
expect a phase transition at the line g = ^(3, below which phase II does not necessarily exist. 
Interestingly, the same linear phase boundary between phases II and III also appeared in the 
numerical results of [18] . The fact that our phase boundary only holds for large \r\ is due 
to our approximation of the exponential of the kinetic term of the original model to second 
order. We expect higher orders to provide the necessary further corrections. Also, higher 
order contributions will yield terms in the potential of the one-dimensional matrix model 
which are of order tr ($"') with n > 4. These terms would give rise to more than two wells 
in the potential leading to additional phases. 

In summary, we conclude that the phase diagram of our theory is a slight deformation of 
that of matrix quantum mechanics, with the addition of a third phase. The phase boundaries 
between phases I and II and phases II and III are given, respectively, by the two curves 

r^-Q^ln^ and , = hs . (5.7) 

10o3 3 
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Appendix 

A. The elliptic integrals 

Throughout the paper, we encountered elliptic integrals, which can be reduced to the follow- 
ing ones: 

Xi= rdA^(62-A2)(A2-a2) , J, = t d\ ^ {h^ - \^){\^ + a^) , 

J a JO 

T2= rdAAV(&'-A2)(A2-a2) , = t d\\^^{b^ - X'){\^ + a^) , (A.l) 

J a Jo 

^3 = I'dX (^(52_A2)(A2_a2))' , Js = d\ (^(52 _ a2)(A2 + a2))' , 

with 6 > a > 0. These integrals can be computed explicitly in terms of elliptic functions. We 
hav^ 

I, = U{{a' + b') Eo - (a2 + 62) - {a' - 52) (F2 - Kq)) , ( A.2a) 

Ji = \{a (-a2 + 62) + a {a^ + 62) Ki) , (A.2b) 

I, = (2 {a' - a'b^ + 6^) Eq - 2 (a^ - a262 + 6^) E2 , ^ , 

15 ' ^ (A. 2c) 

+ (a" - 3a^b^ + 26^) (F2 - Kq)) , 
J2 = —a (2 (a^ + a262 + 6^) Ei - {2a^ + 3a^^ + b") Ki) , (A.2d) 

X3 = If, (2 (a^ - 5a^62 - baH^ + 6^) Eq - 2 (a^ - 50^6^ - 50^6^ + 6^) E2 

35 ^ ^ ' ^ ' (A.2e) 

+ {a^ + 80^6^ - lla26^ + 26^) (F2 - Kq)) , 

J^ = — (2a (-a^ - 50^62 + 50^6^ + 6^) Ei + a [a" + 62) (2a^ + 9a262 - 6^) Ki) , (A.2f) 
35 



^Our conventions for elliptic functions agree with those of Mathematica. 
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where we abbreviated complete and incomplete elliptical functions as follows: 



E, 



•= 



^4 



Kn :=K 



El 



Ki 



-4 



K 



Eo 



E I arcsin ( — 
a 



F2 := F ( arcsin ( — 
V V a 



62 
62 



(A.3) 



B. Further integrals 

The integrals given below are used in the analysis of phase III. Recall our approximation for 
the collective field: 



0o(A) 



-Jpo - a{\ - XminY forAeX, 
vr 



(B.l) 



else . 

The zeros of this collective field define X = [Xl, Xr] and are located at 



— Amin — V — and Ar — Amin + V — 

V a \ a 



(B.2) 



Assuming that the expression under the square root in (f)Q is positive for A G [Al,A/j], we 
have: 



^° ' dA0o(A)(po - a(A - Xrau.?? " 



^"dA0o(A)- „ ^ , 

'"dA<Ao(A)A = ^2^ , dA0o(A)A 



2 _ Po(po + 4aA^i^)) 
8a3/2 



(B.3) 
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